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ABSTRACT 

This paper presents CCD multicolour photometry for the old open cluster NGC 188. The observations were 
carried out as a part of the Beijing-Arizona-Taiwan-Connecticut Multicolour Sky Survey from 1995 February 
to 2008 March, using 15 intermediate-band filters covering 3000-10000 A. By fitting the Padova theoretical 
isochrones to our data, the fundamental parameters of this cluster are derived: an age of f = 7.5 + 0.5 Gyr, a 
distant modulus of (m - M)o = 11.17 + 0.08, and a reddening of E(B -V)- 0.036 + 0.010. The radial surface 
density profile of NGC 188 is obtained by star count. By fitting the King model, the structural parameters of 
NGC 188 are derived: a core radius of Rc - 3.80', a tidal radius of R, = 44.78', and a concentration parameter 
of Co = log(R,/Rc) - 1.07. Fitting the mass function to a power-law function 0(m) oc the slopes of mass 
functions for different spatial regions are derived. We find that NGC 188 presents a slope break in the mass 
function. The break mass is moreak = 0.885 Mq. In the mass range above moreak, the slope of the overall region 
is O' = -0.76. The slope of the core region is a = 1.09, and the slopes of the external regions are a = -0.86 and 
a = -2.15, respectively. In the mass range below mbreak, these slopes are a = 0.12, a = 4.91, a - 1.33, and 
a - -1.09, respectively. The mass segregation in NGC 188 is reflected in the obvious variation of the slopes 
in different spatial regions of this cluster. 

Subject headings: open cluster and association: individual (NGC 188) - stars: fundamental parameters - stars: 
luminosity function, mass function 

1. INTRODUCTION 

Open clusters are important tools to investigate the Galaxy 
evolution and stellar evolution. Due to their accurate param¬ 
eters (e.g. age, metallicity, and distance), open clusters can 
probe th e Galactic d i sk an d trace the Galaxy evolution. For 
example. lChen et ^(l2003h compiled the most complete open 
cluster sample with metallicity, age, and distance data, as well 
as kinematic information, available, and presented some sta¬ 
tistical analysis on spatial and metallicity distributions. Their 
result supports the inside-out Galactic disk formation mech¬ 
anism, in which the invok ed star forma t ion ra te and infall 
timescale vary with radius. iBonatto et al.l (120061) used a sam¬ 
ple consisting of 654 open clusters to determine the thin- 
disk scale height, displacement of the Sun above the Galac¬ 
tic plane, scale length and the open cluster age-distribution 
function. The dynamical state of open clusters is complexly 
determined by the primordial gas expulsion, s tellar evolution , 
dynamics, and the effect of Galactic t idal field (lKroupal2001at 
iWu et al.l 120071: IBonatto et al.ll2005h : (1) the primordial ex¬ 
pulsion of gas unbinds most of cluster, but a significant frac¬ 
tion of it can condense by two-body interactions to became 
an open cluster; (2) stellar evolution at the early phase causes 
large mass loss that will result in the disruption of open clus¬ 
ters; (3) internal dynamical evolution will result in mass seg¬ 
regation and evaporation of low mass star, and (4) the Galactic 
tidal force will strip the stars from the open cluster. Therefore, 
the detailed study of the dynamical state of open clusters can 
examine the stellar evolution theory and Galactic dynamics 
as well as map the Galactic tidal field. The accurate and pre¬ 
cise fundamental parameters, structural parameters and mass 
function are the cornerstone for these investigations. 

The photometric system of the Beijing-Arizona-Taiwan- 
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Connecticut (BATC) Multicolour Survey consists 15 
intermediate-band filters which cover the wavelength range 
3000 - 10000 A. This system is designed to avoid most of 
the known bright and variable night-sky emission lines, and 
it is suitable to derive spectral energy distributions (SEDs) 
of objects. Before February 2006, a Ford Aerospace 2k x 2k 
thick CCD camera was applied, which has a pixel size of 15 
Jim and a field of view of 58' x 58', resulting in a resolution 
of 1".67 pixel '. After February 2006, a new 4k x 4k CCD 
with a pixel si ze of 12 um w as used with a resolution of 
1".36 pixel * (iFan et al.l 120091) . Due to the large field of 
view and deep photometry of the BATC photometric system, 
it is suitable for investigating open clusters in the Milky 
Way. The BATC photometric system has been used to study 
the membership determination, fundamental parameters, 
structural parameters and mass fun ction of open clusters such 
as M67, M48, an d NGC 7789 (iFan etal.l ll9^ IWu et al.l 
I2005L120061 12003). In this paper we will study an old open 
cluster NGC 188 based on the BATC photometric data. 
We will derive the fundamental parameters and structural 
parameters, and discuss the mass function and dynamical 
state for this open cluster. 

The old open cluster NGC 188 has once been regarded 
as the oldest open cluster in the Milky Way. Old open 
clusters are important to understand the early history of the 
Milky Way. Owing to its special location / = 122.85° and 
b - +22.39°, NGC 188 is less contaminated and easy to 
be observed. Therefore, since from ISandagel (119621) many 
works about NGC 188 have been done. Works about the fun- 
damental paramete rs of this cluster are listed in Table 1 of 
iFornal et akl (120071) and Table 1 of this paper. From these ta¬ 
bles we can see that many methods are used to derive the fun¬ 
damental parameters, and the results are somewhat different. 
The age ofNGC 188 is from 4.3 Gyr to 16 Gyr, the metallicity 
is from [Fe/H] = -0.6 to [Fe/H] = 0.69, the reddening value 
is from E{B - V) = 0 to E{B -V) - 0.15, and the distance is 
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from 1275 pc to 2047 pc. We consider that, one hand, these 
scatters result from the reddening-metallicity degeneracy. On 
the other hand, it reflects that different stellar evolutionary 
models and different photometric data will derive different re¬ 
sults. In addition, quite some works adopted solar metallicity 
to NGC 188 and derived the other parameters based on this 
metallicity. Based on the BATC photometric data, we will 
try to find a method to avoid the reddening-metallicity de¬ 
generacy. The color-color diagram based on the photometric 
data in the fo ur BATC bands perf orms well to this purpose 
for NGC 188. iBonatto et al.l (l2005h has derived the structural 
parameters and mass function (MF), and studied the dynami- 
cal evolution of NGC 1 88 i n detail. In addition, t ill now only 
IBonatto & Bical (l2005l) and IBonatto et al.l (l2005h studied the 
MF bre ak of NGC 188, howev er, b ecause of not deep p ho¬ 
tometry. [Bonatto & Bical (l2005h and iBonatto et alJ (l2005l) did 
not observe the MF break of NGC 188. So, in this paper, we 
will study the dynamical evolution of NGC 188 based on our 
deeper photometry. 

We describe our new observational data and photometric re¬ 
duction of NGC 188 in Section 2. In Section 3 we derive the 
fundamental parameters of NGC 188. The structure of NGC 
188 is discussed in Section 4. The MF and the mass segre¬ 
gation are discussed in Section 5. A summary is presented in 
Section 6. 

2. OBSERVATION AND PHOTOMETRIC DATA 

Our observations are based on the BATC photometric sys¬ 
tem. The BATC photometric system utilizes the 60/90 cm 
f/3 Schmidt telescope which is at the Xinglong Station of the 
National Astronomical Observatories, China Academy of Sci¬ 
ences (NAOC). Figure 1 shows the image of NGC 188 in the 
B ATC i band. The deta ils of the BATC sy stem are described 
in lYan et^ (l2000ll and IZhou et alJ (1200 ill . In Table 2, we list 
the corresponding effective wavelength and full width at half 
maximum (FWHM) of each filter, and exposure time and the 
number of frames observed in each filter for our observation. 

Data reductions are carried out by two procedures called 
PIPELINE I and PIPELINE II. These procedures are devel¬ 
oped as standard for the BATC photometric system. Bias 
subtraction and field flatting are carried out by PIPELINE I. 
Instrumental magnitudes of point sources are determined by 
PIPELINE II. After that, the instrume ntal magnitudes ar e cal¬ 
ibrated to the BATC standard system (IZhou et aljl2003h . Eor 
each filter, the average calibration error is below 0.02 mag. 
Eor every image, the mean of the EWHM of the point spread 
function (PSE) is 4.39", which critically samples (« 2.5 pix¬ 
els) the PSE. So, the effect of undersampling on our photo¬ 
metric accuracy is not critical. Eor each star observed more 
than once in a BATC band, the final photometric result in that 
band is the weighted mean. 

3. FUNDAMENTAL PARAMETERS FOR NGC 188 

We ad opt the proper-moti on membership probabilities in 
Table 3 of lPlatais et al.l(l2003h for NGC 188. Matching this ta¬ 
ble and our star catalog, we pick up the common stars, which 
are plotted in Eigure 2. In this paper, stars whose proper- 
motions membership probabilities greater than 60% are used 
to derive the fundamental parameters. We used the photo¬ 
metric data in the seven BATC bands to construct four color- 
magnitude diagrams (CMDs). These combinations are {c- p) 
vs. c, {d - g) vs. d, if - i) vs. /, and (c - m) vs. m. We fit the 
theoretical BATC isochrones (see Section 3.2) to our observa¬ 
tional data of NGC 188 to derive the fundamental parameters. 


3.1. Isochrones of stellar evolutionary models 

Stellar evolutionary models from the Padova group 
(iBertelli et al.lll99^ iGrrardi et alJl200(i I2002L and references 
therein) are widely used. In the Padova stellar evolutionary 
models, iGirardi et al.l (l2002h provided tables of theoretical 
isochrones in such photometric systems as ABmag, STmag, 
VEG Amag, and a standa rd star system. The complete data¬ 
base (IGirardi et al.l 1200^ covers a very large range of stel- 
lar masses (typical l y from 0.6 to 120 Mq). These models 
(IGirardi et al.ll2002l |2008^ are computed with updated opaci¬ 
ties and equations of state, and moderate amount of convec¬ 
tive overshoot. However, the isochrones are presented for 
only 6 initial chemical compositions: [Ee/H] = -2.2490, 
-1.6464, -0.6392, -0.3300, -(-0.0932 (solar metallicity), and 
-1-0. 5595, which are evide ntly not dense enough. It is fortunate 
that iMarigo et alJ (120081) provide tables for any intermediate 
value of age and metallicity via an interactive web interface^. 

3.2. Isochrone fitting 

To determine the main characteristics of the population 
in NGC 188, we fit isochrones to the cluster CMD. We 
used the Padova the oretical isochrones in the BATC system 
(IMarigo et al.ll2008h . Via an interactive web interface (see 
Section 3.1), we can construct a grid of isochrones for dif¬ 
ferent values of age and metallicity, photometric system, and 
dust properties. We use the default models that involve scaled 
solar abundance ratios (i.e., [a/Ee] = 0.0). In pe rforming, 
the C habrier lognormal initial mass function (IME) (IChabrieii 
120011) is adopted. 

3.2.1. Reddening 

To determine the reddening value of an open cluster, a 
method based on the color-color diagram fitting is usually 
used. However, in this method, the degeneracy between red¬ 
dening and metallicity requires that the metallicity is pre¬ 
viously known. Eor instance, adop ting the metallicity of 
[Ee/H] = -0.04 + 0.05 for NGC 188. lSaraiedini et all 
obtained NGC 188’s reddening to be E{B - V) - 0.09 + 
0.02 based on the UBV color-color diagram and the Hyades 
fiduc ial sequence . Ass uming a solar metallicity for NGC 
188, lEornal et~^ (120071) obtained NGC 188’s reddening to 
be E{B - V) - 0.025 + 0.005 based on the (g' - /) and 
(m' - g') color-color diagram and the IGirardi et al.l (120041) 
SDSS isochrones. 

We derive the reddening of NGC 188. Table 1 of 
lEornal et^ (120071) and Table 1 of this paper list the metallic- 
ities of NGC 188 obtained from spectral analysis, which are 
between [Ee/H] = -0.12 and [Ee/H] = h- 0.12. We adopted 
these two values of metallicity to derive the reddening of NGC 
188 on the (d - g) and (/ - i) color-color diagram. In addition, 
the resent results showed that the age of NGC 188 is between 
6.0 and 8.0 Gyr (see Section 3.3 for details). In this paper, the 
reddening E(B - V) is transformed to each BAT C band using 
the extinction coefficient derived b viChenI ( 120001) based o n the 
procedure given in Appendix B of iSchlegel et al.l (119981) . 

Eigure 3 shows the result obtained for three values of red¬ 
dening E(B - V) = 0.0, 0.036, and 0.082 ([Ee/H] = -0.12). 
In the left and right panels, the BATC isochrones are for 6.0 
and 8.0 Gyr, respectively. In Eigure 3, the two vertical dashed 
lines represent the main-sequence band between (d—g) - 0.65 
and 1.05, which is the region where the data and isochrones 

^ http://stev.oapd.inaf.it/cmd 
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Fig. 1. — An image in the field of NGC 188 observed by the BATC i band. The image size is 58' X 58'. 


should fit best. From Figure 3, we can see that the best fit 
of the BATC isochrone to the data points was obtained for 
E(B -V)- 0.036. Figure 4 shows the result as Figure 3 did, 
but for [Fe/H] = +0.12. From Figure 4, we can also see that 
the best fit of the BATC isochrone to the data points was ob¬ 
tained for E{B - V) = 0.036. As a result, we arrived at a 
reddening value of E{B - V) - 0.036 + 0.010 for NGC 188, 
where the quoted uncertainty only includes an estimate of the 
error in the reddening determination. 

3.2.2. Metallicity, age, and distance modulus 

Since the value of reddening has been fixed, we go ahead 
to determine the other parameters of NGC 188 in the CMDs. 
Because of the age-metallicity degeneracy, we cannot simul¬ 
taneously determine age and metallicity by the fitti ng method. 
We ad opt a value of metallicity [Fe/H] = 0.12 in iFriel et alJ 
(120101) . which was obtained from spectral analysis of red gi¬ 
ants of NGC 188. Now, fitting the theoretical isochrone to our 
color-magnitude data points became straightforward. First, 
we chose the age for which the fitted isochrone had the same 
color of the turn-off as our data points. Then, by shifting the 
isochrone up and down (changing the distance modulus), we 
finally obtained the best fit. 

Based on the age obtained by the previous works (see Ta¬ 
ble 1 of this paper and Table 1 of Fornal et al. 2007 in de¬ 
tail), we used the interactive web in Section 3.1 to construct 
a fine grid of isochrones about ages, sampling an age range 
6.0 < t < 8.0 Gyr at intervals of 0.1 Gyr. The metallic¬ 
ity [Fe/H] = log(Z/Zo) = 0.12 where Zq ~ 0.0152, so this 
metallicity corresponds to Z = 0.0200. Using this fine grid of 
isochrones, we obtained the fundamental parameters of NGC 
188. A set of the best fitting parameters for our data is adopted 
as follows: an age of f = 7.5 + 0.5 Gyr, and a distance modu¬ 
lus of im - M)q - 11.17 + 0.08. The error of the age of NGC 
188 is estimated by eye in the CMDs. Figure 5 shows a set of 
CMDs in the BATC system for NGC 188. From Figure 5, we 
can give an estimate of the error of age to be ±0.5 Gyr. The 


error of the distance modulus is also estimated by eye in the 
CMDs. Figure 6 shows a set of CMDs in the BATC system 
for NGC 188. From Figure 6, we can give an estimate of the 
error of distance modulus to be ±0.08. 

3.3. Comparing with previous works 

We compare our best fitting parameters for NGC 188 with 
those of previous studies in the past two decade. For red¬ 
dening E(B - V), many works (see Table 1 of Fornal et al. 
2007) obtained the values of reddening to be larger than 0.08. 
Our best fitti ng result of 0.036 ± 0 .010 is in very agreement 
with those in ICarraro et al.l (Il994h . who obtained the values 
of reddening to be 0.03 and 0.04. Our result is also consis¬ 
tent with that of iFornal et al.l ([2003), who obtained the value 
of reddening to be 0.025 ± 0.005. For age, o ur best fitting 
of 7.5 ± 0.5 Gyr is consisten t wi th those of ICarraro et aH 
(ll994ft . [Sarajedini et al.l (I1999^ . and iFornal et al.l(120071) . who 
obtained the ages of NGC 188 to be 7.5 Gyr and 7.0 Gyr, 
7.0 ± 0.5 Gyr, and 7.5 ± 0.7 Gyr, respectively. For dis¬ 
tance modulus (m - M)q, our best fitting of 11.17 ± 0.08 
or 17 1 4 ± 64 pc) is consiste nt wi th those of ICarraro et al.] 
1994h . [^aiedini et al.l(ll999h . and iFornal et al.l (l2007h . who 
obtained the distances ofNGC 188 to be 1705 pc and 1720 pc, 
1710 ± 70 pc, and 1700 ± 100 pc, respectively. 


4. STRUCTURAL PARAMETERS OE NGC 188 

The structure of a cluster is described by means of the ra¬ 
dial surface density profile, which is defined as the projected 
number of stars per area in the direction of the c luster. In this 
paper we adopt the empirical formula derived bv iKingl (119621) 
to fit the radial surface density profile of NGC 188. 

In order to enlarge the number of member stars of NGC 
188, w e used the SED fitting method developed by IWu et’aH 
(l2006l) to search more member stars. For the jth star, a 
parameter S can be defined: 
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Fig. 2. — Observed (c - p) vs. c, {d - g) vs. d, (/ — i) vs. /, and {c — m) vs. m CMDs of NGC 188. Dots represent the stars with the membership probability 
of 0% < P < 60%, crosses represent the stars with the membership probability of 60% < P < 90%, and diamonds represent the stars with the membership 
probability of 90% < P < 100%. 


Sj{t,Z,(m-M)o,E(B-V)) 


= Z 


[mu - Mi(t, Z, (m - M)o, E(B - V))Y 




( 1 ) 


where M,(f, Z, (m-M)o, E{B-V)) is the theoretical magnitude 
in the /th BATC band, corrected by the distance modulus {m - 
M)o and reddening E{B - V) and computed from the chosen 
theoretical isochrone model with the age t and metallicity Z. 
The values of (m - M)q, E(B- V), t, and Z have been obtained 
in Section 3. Here mij and cr,y are the observed magnitude 
and its error, respectively, of the y'th star in the /th band, and 
n is the total number of observed bands for the y'th star. For 
Mi with different stellar masses, the minimum of S j, S 
can be obtained for the y'th star with the chosen theoretical 
model. If the observed SEDs can match the theoretical SEDs, 
the parameter S min should be the distribution with n - P 
degrees of freedom, where P is the number of free parameters 
to be solved. The integral probability at least as large as S p„i„ 
in the x^ distribution with n - P degrees of freedom is taken 
as the “photom etric" membership probability of the y'th star 
(IWu et alj|2006ll . We determine 2131 stars to be photometric 
members with Pphot greater than 50%. The d-g versus d 
CMD of these stars are shown in Eigure 7. 

We used these 2131 member stars to determine the radial 
surface density profile of NGC 188. We adopted the cen¬ 
tral position for NGC 188 at a n - 00^47™53'^ and do - 
-i-85°15'30" (J2000.0) following lonitoetin (Hop). The 
radial surface density profile is determined by counting stars 
inside concentric annuli with a step of 1.0' in radius of 0' < 


r < 20' and counting stars inside concentric annuli with a 
step of 2.0' in radius of 20' < r < 30'. The radial surface 
density profile of NGC 188 is shown in Eigure 8. The fit 
has been performed using a non-linear least-square fit routine 
which uses the 1 cr Poisson errors as weights. The fitted struc¬ 
tural parameters obtained here are as follows: a core radius of 
Rc - 3.80', a tidal radius of R, - 44.78', and a concentration 
parameter of Co = log {RtIRc) = 107. With the distance to 
the Sun d^ - 1714 pc obtained in Section 3, the core radius 
and tidal radius of NGC 188 turn out to be = 1 -89 pc and 
Rf = 22.33 pc, respectively. The fitting curves are shown in 
Eigure 8 with solid lines. As shown in Eigure 8, the radial 
surface density profile is reproduced well by the King model 
within uncertainties. 

We compare our structural parameters for NGC 188 with 
those of previous studies. Eor the core radius Rc, our 
best fitting result of 3.80' is somewhat larger than that of 
iBonatto et al.l (l2005h . who obtained the value of core radius 
to be 3 .1 +0.2'. However, ou r result is somewhat smaller than 
that of IChumak et al.l (120101) . who obtained the value of core 
radius to be 4.45'. The tidal radiu s of R, - 44.78' ob tained 
here is somewhat larger than that in iKeenan et al.l (119731) . who 
obtained the tidal radius of N GC 188 to be 27.0 + 5.0'. Our 
result is consistent with that in IBonatto et al.l (1200.51) . who ob¬ 
tained the tidal radius of NGC 188to be44:t :9'. And our result 
is smaller than that in IChumak et al.l (1201 Ol) . who derived the 
tidal radius of NGC 188 to be 71.6'. 


5. MASS FUNCTION OF NGC 188 
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Fig. 3. — Observed d - g vs. f — i color-color diagram for the main-sequence stars with membership probability larger than 60% in NGC 188. The dotted, solid, 
and dashed lines are BATC isochrones ([Fe/H] = -0.12) with three different reddening values E{B -V) = 0.0, 0.036, and 0.082. In the left and right panels, the 
BATC isochrones are for 6.0 and 8.0 Gyr, respectively. The two vertical dashed lines represent the main-sequence band between (d - g) = 0.65 and 1.05, which 
is the region where the data and isochrones should fit best. Error bars represent the photometric en'ors. 



Fig. 4. — Observed d - g vs. f — i color-color diagram for the main-sequence stars with membership probability larger than 60% in NGC 188. The dotted, solid, 
and dashed lines are BATC isochrones ([Fe/H] = +0.12) with three different reddening values E{B — F) = 0.0, 0.036, and 0.082. In the left and right panels, the 
BATC isochrones are for 6.0 and 8.0 Gyr, respectively. The two vertical dashed lines represent the main-sequence band between {d - g) = 0.65 and 1.05, which 
is the region where the data and isochrones should fit best. EiTor bars represent the photometric errors. 


Since a large field of view of NGC 188 is observed here, we 
have an opportunity to study the MF and the mass segregation 
for NGC 188. The member stars of NGC 188 obtained here 
by the SED fitting method is used. In addition, we take the 
stars below the main-sequence turn-off to derive MF. We used 
stars with Pphot greater than 60% to avoid the field-star con¬ 
tamination. In order to study the mass segregation of NGC 
188, we must know the mass functions in different regions. 
We select the core region (O' < R < 3.80'), 3.80' < R < 15', 
15' < R < 30', and the whole region (O' < R < 30'). Com¬ 
pleteness estimation is necessary to obtain the MF of NGC 
188. We used artificial-sta r tests on our da ta to determine 
completeness corrections as iWu et alJ (I2007h did. The com¬ 
pleteness correction for the whole region and for three differ¬ 
ent regions of NGC 188 are listed in Table 3 for the BATC e 
band. In the mass functions we included only the values for 
which the completeness corrections are 0.5 or higher for all 
regions in NGC 188. Therefore, limiting magnitude e - 20.0 
is adopted. After correcting for completeness, the mass func¬ 
tions for different regions of NGC 188 were obtained, which 


are shown in Figure 9. 

In Figure 9, the MF break is obvious. The break mass is de¬ 
termined to be near mbreak = 0.885 Mq. The MF break of clus¬ 
ters essentially reflects the internal dynamics of clusters on 
the mass function and/or some fundamental properties of the 
initial mas s function associated to different conditions in star 
formation (iBonatto & Bicall2005h . The power-law function 
(f>(m) oc is used to fit the MF. The fit has been performed 
using a non-linear least-square fit routine which uses the 1 cr 
Poisson errors as weights. The best fitting power-law func¬ 
tions are shown in Figure 9 with solid lines, and the slopes 
are listed in Table 4. In the fitting of the core region, we did 
not use the point which only includes one member star. The 
slope of MF is fiat in the core region and become steeper as 
the distance increases, which reflects the mass segregation of 
NGC 188. The overall MF slope a = -0.76 in the mass range 
0.885 Mq < M < 1.125 Mq is much flatter than the universal 
IMF a - -2.2 + 0.3 in the mass range 0.5 Mq < M < 1.0 Mq 
(lKrour)all2001bl I2002h . This flat slope reflects the advanced 
dynamical state of NGC 188. The number of stars for NGC 
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Fig. 5.— Observed {c - p) vs. c, (d — g) vs. d, (/ — i) vs. /, and {c - m) vs. m CMDs of NGC 188 for stars with membership probability larger than 60%. The 
blue, red, and green lines are BATC isochrones (Z = 0.0200, (m - M)q = 11.17, and E{B — V) = 0.036) with three different ages t = 7.0 Gyr, 7.5 Gyr, and 8.0 
Gyr. En'or bars represent the photometric errors. 


188 can be estimated by the power-law function. We derived 
the number of stars in NGC 188 based on the power-law func¬ 
tion obtained here to be 5776. Based on the number of stars 
obtained here, we can estimate the relaxation time of NGC 
188. The relaxation time is calculated according to this for¬ 
mula freiax = A^fcross/CSltiV), where fcross = R/o'v is the cross 
time, N is the number of stars, and cr,, is the velocity disper¬ 
sion (iBinnev & Tremai n3l2008h. We ado pted a typical cr„ s; 
3 kms~^ from iBinnev & Merrifieldl (Il998h . and the relaxation 
time of NGC 188 is freiax = 595.6 Myr. The ratio of the cluster 
age to the relaxation time is t = Age/freiax = 12.6. Therefore, 
NGC 188 is in the advanced dynamical state. Many number 
of low mass stars have been evaporated and stripped from this 
open cluster. 

6. SUMMARY 

In this paper, we present photometry of NGC 188 in 15 
intermediate-band filters of the BATC system. By fitting the 
theoretical isochrones to the observational data points, we de¬ 
termined three fundamental parameters for NGC 188: a red¬ 
dening of E{B -V)- 0.036 + 0.010, an age of t - 1.5 ± 0.5 
Gyr, and a distance modulus of (m - M)q - 11.17 + 0.08. 

Bas ed on the SED fitting method developed by IWu et alJ 
daooS), we obtained a sample of NGC 188 member stars con¬ 
sisting of 2131 stars with Pphot greater than 50%. We deter¬ 
mined the radial surface density profile of NGC 188 using this 


sample. By fitting the empirical density law of lKind (Il962h 
to our radial surface density profile, the structural parame¬ 
ters of NGC 188 are derived: a core radius of Rc - 3.80', a 
tidal radius of Rt - 44.78', and a concentration parameter of 
Co = 1.07. 

We studied the details of mass function of NGC 188 based 
on the photometry which are obtained in a large field of view 
and deep observation. The ME is fitted with a power-law func¬ 
tion, and the slopes of mass functions in different regions are 
derived. The ME slopes flatten from the outskirts to the in¬ 
ner regions which indicates mass segregation in NGC 188. 
The ME break in NGC 188 is observed, and the break mass is 
wibreak = 0.885 Mq. We also derived the number of stars and 
the relaxation time. The mass segregation and the large ratio 
of the cluster age to the relaxation time indicate that NGC 188 
has gone strong dynamical evolution. 
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Fig. 7.— Observed d~gvsd CMD of NGC 188 for member stars which are determined by the SED htting method. Error bars represent the photometric errors. 
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Fig. 8.— Radial surface density profile of NGC 188. The top panel uses logarithm coordinate, and the bottom panel uses linear coordinate. Error bars represent 
the 1 cr Poisson errors. The fitting parameters of King model are adopted as follows: a core radius of Rc = 3.80' and a tidal radius of R[ = 44.78'. 
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Fig. 9.— Mass functions in different spatial regions of NGC 188. Error bars represent the 1 cr Poisson errors. Solid lines are the best fitting power-law functions. 
The spatial range and the slope of power-law function are given in each panel. We used stars with Pphot greater than 60% to avoid the field-contamination. 
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TABLE 1 

Literature estimates of fundamental parameters for NGC 188 


Reference 

Distance Modulus 

Distance 

(pc) 

Age 

(Gyr) 

E{B - V) 
(mag) 

[Fe/H] 

Technique 

Fornal et al.(2007) 

11.23 ±0.14 

1700 ± 100 

7.5 ± 0.7 

0.025 ± 0.005 

O.OO(adopted) 

u g r i' z' photometry,Girardi isochrone fits 

Meibom et al.(2009) 

11.24 ±0.09 

1770 ± 75 

6.2 ± 0.2 

0.087 (adopted) 

-0.1 - O.O(adopted) 

Analysis of eclipsing binaries 

Friel et al.(2010) 





0.12 ±0.02 

High-resolution spectral analysis of red giants 

Jacobson et al.(2011) 





-0.03 ± 0.04 

High-resolution spectral analysis 

This work 

11.21 ±0.08 

1714 ± 64 

7.5 ±0.5 

0.036 ±0.010 

0.12(adopted) 

BATC multi-color photometry, isochrone fits 
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TABLE 2 

Parameters of 15 BATC filters and statistics of our observations 


Number 

Filter Name 

/leff 

(A) 

FWHM 

(A) 

Exposure 

(s) 

Number of Images 

1 

a 

3360 

360 

4800 

7 

2 

b 

3890 

340 

11010 

18 

3 

c 

4194 

309 

7020 

8 

4 

d 

4540 

332 

6840 

7 

5 

e 

4925 

374 

23940 

24 

6 

f 

5267 

344 

11709 

28 

7 

g 

5790 

289 

7620 

14 

8 

h 

6074 

308 

6780 

9 

9 

i 

6656 

491 

5800 

16 

10 

j 

7057 

238 

7080 

10 

11 

k 

7546 

192 

9300 

11 

12 

m 

8023 

255 

12780 

17 

13 

n 

8484 

167 

12780 

10 

14 

o 

9182 

247 

17985 

15 

15 

p 

9739 

275 

19800 

19 


TABLE 3 

Completeness analysis results for NGC 188 


Ae 

0.0' < r < 30' 

0.0' <r< 3.80' 

3.80' < r < 15 

15 < r < 30 

13 

- 14 

0.91 

0.37 

0.92 

0.91 

14 

- 15 

0.92 

0.56 

0.93 

0.92 

15 

- 16 

0.96 

0.96 

0.97 

0.97 

16 

- 17 

0.91 

0.90 

0.92 

0.92 

17 

- 18 

0.83 

0.82 

0.84 

0.84 

18 

- 19 

0.68 

0.69 

0.69 

0.69 

19 

-20 

0.53 

0.64 

0.56 

0.56 

20 

-21 

0.13 

0.47 

0.14 

0.14 


TABLE 4 

Fitted parameters for mass functions of NGC 188 


Distance r 
(arcmin) 

a (m< Hlbrcak) 

0.6 - 0.885 

a(m> mbreak) 

0.885 - 1.125 

0.0-30.0 

0.12 

-0.76 

0.0-3.80 

4.91 

1.09 

3.80-15.0 

1.33 

-0.86 

15.0 - 30.0 

-1.09 

-2.15 



